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Abstract: Fast scramblers are quantum systems which thermalize in a time scale 
logarithmic in the number of degrees of freedom of the system. Non-locality has been 
argued to be an essential feature of fast scramblers. We provide evidence in support 
of the crucial role of non-locality in such systems by considering the approach to ther- 
malization in a (strongly-coupled) high temperature non-commutative gauge theory. 
We show that non-locality inherent to non- commutative gauge theories does indeed 
accelerate the rate of dissipation in the heat bath in stark contrast to the slow random 
walk diffusive behavior prevalent in local field theories. 
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1. Introduction 

It has been known since the days of the "membrane paradigm" [1] that the event 
horizon of black holes, or more precisely the stretched horizon, is endowed with physical 
properties from the vantage point of some fiducial observers. This is in contrast with 
the free infalling observers which do not experience the horizon. In particular, from 
the viewpoint of these fiducial observers, a charge falling into the horizon is seen as 
having its charge smeared over the horizon in a time that grows logarithmically with 
the Schwarzschild radius of the black hole. 

In a remarkable paper, Preskill and Hayden [2] posed a critical challenge to the 
resolution of the information paradox. In a nutshell, their observation lead to the 
conclusion that the time scale involved in possibly sharing information between two 
observes, one hovering outside the black hole horizon and the other falling through 
the horizon while carrying information, is dramatically shorter than what was believed 
possible until then. The fascinating answer to this challenge, the so-called "fast scram- 
bling" property of the horizon, was conjectured by Sekino and Susskind in [3]. This 
property is tantamount to the extremely rapid spreading of the charge as we alluded 
to earlier. The time scale involved in scrambling grows only logarithmically with the 
entropy. This is a much shorter time scale as a function of the entropy than the time 
that a particle subject to Brownian motion would take to travel through a distance 
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commensurate with the "size" of the system [3]. The time scale in the latter case grows 
with a positive power of the entropy. 

It was argued in [3] that non-locality has to be one of the ingredients that generate 
fast scrambling. Local interactions, on the other hand, would lead to slow diffusion as 
in Brownian motion. Another possibly necessary ingredient to generate the fast scram- 
bling property of the stretched horizon is the presence of strong coupling interactions 
among the horizon's degrees of freedom [3]. 

In this paper, we will emulate non-locality and strong coupling by considering a 
strongly-coupled non-commutative gauge theory at high temperature: four-dimensional 
non-commutative (supersymmetric) Yang-Mills theory, to be precise. As the non- 
commutative system under consideration is strongly coupled, we use the AdS/CFT 
correspondence as a toolbox for analyzing various thermal aspects of the system. We 
show that the rate of decay of a disturbance propagating in this thermal medium is 
parametrically much larger than in the local setting of conventional field theories. This 
will inevitably lead to rapid thermalization 1 . It is interesting to note that computations 
at weak coupling do not exhibit an enhancement of the thermalization rate. This could 
then imply that fast thermalization in non-commutative field theories is a consequence 
of having both strong coupling and non-commutativity. See [4, 5, 6, 7] for other works 
in analyzing the fast scrambling conjecture. 

The existence of a connection between the ultraviolet and the infrared in the non- 
commutative setting has been well studied over the years [8] . One of its manifestation 
is the large transverse size of any degree of freedom moving with a large momentum 
k, where transversality here is with respect to the direction of k. The transverse size 
grows with k when the momentum in the non-commutative direction is larger than 
where 9 is the Moyal area. Indeed, the transverse size scales like Ok. This property 
is at the heart of the rapid decay of high momenta modes into excitations of the heat 
bath. 

The paper is organized as follows: in the following section we briefly review the 
gravity dual of a strongly coupled large- iV non-commutative gauge theory at finite tem- 
perature and then set-up the study of the approach to thermalization in this context. 
In section three, we exhibit the computation of the decay rates in the non-commutative 
heat bath and compare it to the results for the local case of strongly-coupled commu- 
tative gauge theories. More specifically, we compute the decay rates in three cases: 
once when the non-commutative thermal bath is perturbed by an operator dual to a 

x It should be noted that the rapid thermalization that we observe here is for the boundary non- 
commutative thermal bath and does not seem to be directly connected to the properties of the stretched 
horizon of the dual black brane background as observed by a fiducial observer hovering outside the 
horizon. 
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minimally-coupled massless scalar field in the bulk, then with the operator dual to the 
dilaton fluctuation of the bulk, and finally with the operator dual to the bulk axion fluc- 
tuation. Furthermore, we compute the decay rates in the regime of parameters where 
the temperature of the bath is taken to be larger than #~ 1 / 2 along with the momentum 
of the modes in the non-commutative plane being greater than the temperature (which 
necessarily implies that the momentum is also greater than 6~ l l 2 ). This is the regime 
where one expects the deviations from the results of the commutative setting to be 
more pronounced 2 . We close with some concluding remarks and open questions at the 
end. 



2. A Gravity Dual For Non-commutative SYM 

In string theory, non-commutative gauge theories can be easily realized on the world- 
volume of D-branes with a constant Neveu-Schwarz B^ u field provided that one takes 
a special limit to decouple the open and closed string modes [9, 10, 11]. Basically, one 
scales the string tension to infinity and the closed string metric to zero while keeping 
the background B^ u field fixed [11]. For much of the discussion in this paper, we focus 
on the four-dimensional non-commutative Yang-Mills theory at finite temperature. We 
also take the non-commutativity parameter to be non-zero only in the (x 2 , x 3 )-plane, 
that is [x 2 ,x 3 ] ~ id. In the spirit of the AdS/CFT correspondence [12], the dynamics 
of four-dimensional non-commutative Yang-Mills theory at large iV and at strong 't 
Hooft coupling should be described by a bulk gravitational system. Indeed, a proposal 
for the dual gravity background was given in [13, 14] which, in the string frame, reads 



ds 2 = L 2 



2 



u 2 f(u)dt 2 + u 2 dx\ + u 2 h(u) (dx\ + dxf) + + dttl 



B 23 = L 2 a 2 u 4 h{u), e 2 * = g 2 h(u), (2.1) 

Coi = L 2 a 2 g- X u\ F 0l23u = AL^g^Ku), 



where L 4 = AngNa' 2 with g being the value of the string coupling in the IR, and 



1 7/ 4 

h(u) = —, f(u) = 1 - ^, (2.2) 

V ; l + a 4 w 4 ' JK J u K J 

where Uh is the horizon radius. 

2 This is due to the fact that at temperatures larger than the typical momenta in the heat 

bath is of the order of temperature and hence the transverse size of the modes of the heat bath is 
large and gets larger as the temperature increases. One therefore expects the non-local effects to be 
important at those high temperatures. 
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The parameter a which appears in the above expressions is related to the non- 
commutativity parameter 6 of the boundary theory through a = X^-^y/O. Here, A = 
L A /a' 2 is the 't Hooft coupling of the boundary large- N non-commutative Yang-Mills 
theory. For u <C a -1 , the background (2.1) goes over to the AdS 5 x S 5 solution. As 
expected, this observation just reflects the fact that the non- commutative boundary 
theory goes over to the ordinary commutative Yang-Mills theory at length scales much 
greater than X^^yfO. For u ^> a -1 , on the other hand, the background (2.1) shows 
significant deviation form the AdSs x S 5 solution. In particular, the bulk spacetime 
is no longer asymptotically AdSs. The boundary theory interpretation of this regime 
just means that the effect of non-commutativity becomes pronounced for length scales 
which are at the order of, or smaller than, A 1 / 4 ^ 3 - To summarize, the background 
(2.1) represents a flow in the boundary theory from a UV fixed point which is non- 
commutative Yang-Mills at large N and large A to an IR fixed point given by the 
ordinary commutative Yang-Mills (again, at large iV and large A). 

Also, the Hawking temperature of the above solution, which is interpreted as the 
temperature of the non-commutative boundary theory, is given by T = u^/ir. Notice 
that this temperature is the same as the temperature one obtains for the Schwarzschild 
AdSs solution, which is dual to a thermal state of four-dimensional commutative Yang- 
Mills theory. Indeed, it is easy to show that all the thermodynamic quantities obtained 
from (2.1) are the same as the ones obtained from the Schwarzschild AdS 5 solution 
[14]. 

3. Approach to Thermalization 

In the context of the AdS/CFT correspondence [12, 15, 16], turning on the fluctuations 
of the bulk background translates into deforming the boundary field theory by the 
operators dual to the bulk fluctuations. In the case where the background is a large 
black hole, which describes a thermal state in the boundary field theory [17], the bulk 
fluctuations correspond to the perturbations of the thermal bath in the boundary field 
theory. Due to the existence of the horizon in the bulk, the fluctuations of the black hole 
background fall into the horizon. From the perspective of the boundary field theory, 
the decay of the bulk fluctuations into the horizon translates into the decay of the 
perturbations of the thermal bath and their eventual thermalization, after some time 
scale, with the rest of the bath [18]. In other words, the decay of the bulk fluctuations 

3 Since the reliability of the solution (2.1) requires L 4 /a' 2 = A to be large, the effect of non- 
commutativity in the boundary theory is visible even at large length scales [14]. This should be 
compared with the weak-coupling regime of the theory where the threshold length scale (beyond 
which the theory becomes effectively commutative) is roughly at the order of \/6. 
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into the horizon holographically describes the approach to thermal equilibrium in the 
boundary field theory. 

In the linear response theory, the time scale required for a perturbation to ther- 
malize with the rest of the bath can be read off from the imaginary part of the poles 
in the retarded Green's function of the operator by which the thermal bath has been 
perturbed. Indeed, the late-time decay is dominated by the imaginary part of the 
lowest of such poles (in the complex frequency plane). While the computation of the 
thermalization time scales is relatively easy in the perturbation theory [19, 20], it is 
a formidable task when the thermal bath is strongly coupled. Fortunately, this is the 
regime where the AdS/CFT correspondence is most useful. In fact, the poles in the 
retarded Green's function of the boundary theory operators are mapped, through the 
correspondence, to the quasi-normal frequencies of the bulk fluctuations which are dual 
to those operators, a connection which was first proposed in [21] and later sharpened 
in [22, 23]. See [24] for a relatively recent review on the quasi-normal modes of var- 
ious backgrounds (with different asymptotics) from the perspectives of both general 
relativity and the AdS/CFT correspondence. 

The quasi-normal modes of a perturbation of a gravitational background are de- 
fined as the solutions to the equation of motion of that perturbation which are infalling 
near the horizon and normalizable at the boundary (when the background is asymp- 
totically AdS). The quasi-normal frequencies are then the set of complex frequencies 
for which such solutions exist. Thus, for strong-coupling thermal baths, the thermal- 
ization time scales can simply be obtained by computing the quasi-normal modes of 
the dual black hole backgrounds. While there are semi-analytic methods for computing 
special values of the quasi-normal frequencies, their generic values, however, cannot be 
computed analytically for most backgrounds. In such cases, one can compute them 
numerically with high precision using a variety of techniques developed over the years 
in the literature [24]. 

Our goal is to study the approach to thermalization in a strongly coupled non- 
commutative thermal bath, which is dual to the non-extremal background given in (2.1). 
In particular, we would like to analyze to what extent the thermalization properties of 
a non-commutative thermal bath deviate from its commutative counterpart. To do so, 
we will consider three examples of fluctuations of the dual supergravity background: a 
minimally-coupled scalar field, fluctuation of the dilaton as well as the fluctuation of 
the Ramond-Ramond (axion) scalar field. 

3.1 Minimally-Coupled Scalar 

We start by considering perhaps the simplest possibility. That is, we perturb the non- 
extremal background (2.1) by a minimally-coupled massless scalar ip, which we assume 
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to be dual to some scalar operator O in the non-commutative boundary theory. Our 
goal is to calculate the time scale for this perturbation to thermalize with the non- 
commutative bath. As we alluded to the above, we need to compute the quasi-normal 
modes of the minimally-coupled scalar field fluctuations around the background (2.1). 
The equation of motion for (p (with no dependence on the coordinates of S 5 ) reads 



V-9 

which, in the non-extremal background (2.1), takes the form 

= d u [u 5 /(«)M +u[- j^d 2 + c£ + -L + d l 3 ) ] V9 . (3.2) 

We emphasize that the minimally-coupled massless scalar that we are studying 
here is not related to the metric fluctuations of the background (2.1) such as 5g l x i ■ As 
shown in [14], the equation of motion for 5g f xl (with no dependence on S 5 ) satifies the 
same equation as in (3.1) provided that Sg^ depends neither on t nor on x 1 directions. 
It is not clear to us whether the equation of motion for any generic fluctuation of the 
background (2.1), or their gauge invariant combinations, could be given by the equation 
(3.2). As such, the analysis in this section should be regarded as a warm-up exercise. 
Nevertheless, we believe that the lesson we learn by analyzing the case of a minimally- 
coupled massless scalar will prove useful when we consider, in the subsequent section, 
some genuine fluctuations of the background such as the dilaton and axion fluctuations. 

To solve the equation (3.2), we go to momentum space by Fourier transforming 
ip(u, t, x): 

(p(u, t, x) ~ e- iu)t+ll ^{u; u, k). (3.3) 

To avoid clutter in the notation, we have denoted the Fourier modes of the scalar field 
also by tp. The distinction should be clear from the arguments of ip, or the context 
in which it is being used. In order to see any non-trivial effect on the behavior of the 
modes tp(u; u, k) due to the existence of a non-zero Neveu-Schwarz field in the 
bulk, especially in the asymptotic u — > oo region where the behavior of ip(u — > oo; u, k) 
holographically encodes information about the source and the condensate of the dual 
operator 0(u,k), one should consider the modes with a non- vanishing momentum in 
the (x 2 ,x 3 ) plane. Without loss of generality, we take k 2 = and k 3 ^ and, for 
simplicity, we also consider the modes with k\ = 0. The equation of motion (3.2) then 
becomes 

u 2 k 2 



d u [u 5 f(u)d u ip] + 



u 



/(«) Ku) 



V = 0, (3.4) 
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where k = k 3 . 

Some comments on the boundary theory operator dual to ip(u,t,x) are in order 
here. Normally, in the AdS/CFT correspondence, one associates a gauge invariant 
operator 0(t,x) in the boundary theory to the bulk fluctuation ip(u,t,x) in the sense 
that near the boundary of the spacetime the non-normalizable mode of ip(u — > oo,t,x) 
sources the dual operator 0(t, x) while its normalizable mode gives the vacuum expec- 
tation value of 0(t, x). The situation is less clear when the boundary gauge theory 
is non- commutative. This is partly due to the fact that in non-commutative gauge 
theories, non-commutativity of the space mixes with the gauge transformations and, 
as such, there are no gauge invariant operators in position space. In momentum space, 
however, one can construct gauge invariant operators Oik 11 ) by smearing the gauge co- 
variant operators 0(x^) transforming in the adjoint over an open Wilson line W(x, C) 
of definite size according to [25, 27, 26] 

0(k) = J d 4 x 0(x) * W(x, C) * e ihx , (3.5) 

where * denotes the Moyal product. The expression in (3.5) reduces to the corre- 
sponding ordinary gauge invariant local operators in the IR. In the UV, on the other 
hand, this very simple modification, has a drastic consequence for the behavior of the 
correlation functions of the gauge invariant operators 0{k). 

In the following we assume that the bulk fluctuation ip(u;u,k) is dual to a gauge 

— # 

invariant operator O(to, k) of the form (3.5) in the sense that in the boundary non- 
commutative gauge theory there exists a coupling of the following form in the momen- 
tum space, 

J dcud 3 kip (-to,-k)0(uj,k), (3.6) 

where f (u,k) is the source term which is determined from the asymptotic non- 
normalizable behavior of <p(u,u,k) given some suitable boundary condition near the 
horizon. 

To obtain the retarded Green's function of the dual operator 0(u, k), equation (3.4) 
should then be solved with infalling boundary condition at the horizon [22]. Close to 
the horizon, the solution behaves as 

<p(u) ~ f{u)^/^ T \ u^u h , (3.7) 

where the + and — signs represent outgoing and infalling waves, respectively. So, we 
discard the solution with the + sign. Near the boundary u — A — > oo, the equation 



-7- 



(3.4) can be conveniently cast into a Mathieu differential equation. The asymptotic 
behavior of the solution then takes the form 

<p(u) ~ A(u, k) u~ 5/2 e ka2u + B(co, k) vT^e'^ . (3.8) 

The first term in the above equation blows up as the cutoff A — > oo. So, in the spirit of 
the AdS/CFT correspondence, the source for the operator dual to (f in the boundary 
theory should be read off from that term. However, extracting the source of the dual 
operator from (3.8) is ambiguous. The ambiguity arises because, as seen from (3.8), 
for any cutoff A the relation between the would-be source A(u, k) and the asymptotic 
behavior of tp involves a function of k, namely e . As a result, the normalization 
of the dual operator, and hence the normalization of the correlation functions, in the 
boundary theory is not fixed 4 . This feature, which also arises in [28, 29], can ultimately 
be attributed to the non-local nature of the boundary theory in the UV. As advocated 
in [14] (see also [30]), a reasonable way forward is to define the source for the dual 
operator from the asymptotic behavior of the bulk fluctuation in such a way that the 
correlation functions of the operator reproduce the ordinary commutative results in the 
IR. Nevertheless, regardless of how the normalization of the boundary theory operator 
is fixed, the ratios of the correlation functions are unambiguous [14, 27]. Also, note 
that the overall /c-dependent normalization should not affect the poles (in the complex 
frequency plane) of the retarded Green's function of the dual operator. 

In order to find the poles in the retarded Green's function of the dual operator 
0(u, k), we numerically solve (3.4) with infalling boundary condition near the horizon 
and normalizable boundary condition near the boundary 5 . In other words, we compute 
the quasi-normal frequencies of f in the non-extremal background (2.1). To proceed, 
we find it convenient to define a new radial coordinate z = Uh/u. In terms of this new 
radial coordinate, the horizon is now at z — 1 while z — > is the asymptotic boundary 
of the spacetime. Using the relation Uh = 7rT with T being the Hawking temperature 
of the background, the equation of motion (3.4), expressed in terms of the new radial 

4 This is in contrast to the usual examples in the AdS/CFT correspondence where the background 
is asymptotically AdS and the relation between the asymptotic behavior of a bulk fluctuation and the 
source of the corresponding dual operator does not depend on k. Indeed, such a relation is fixed by 
the conformal symmetry of the boundary theory. 

5 The normalizable boundary condition near the boundary seems to be a natural choice given that 
the potential term in equation (3.4), once written in the form of a Schrodingcr equation, blows up near 
the boundary, see the equation (5.2) in the appendix. This is similar to the case where the background 
is asymptotically AdS where the potential term in the equation of motion for a minimally-coupled 
scalar also blows up at the boundary of the spacetime. 
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coordinate, takes the form 

<p(z), (3.9) 



= f(z)zV(z) + [zf'(z) - 3f(z)} z<p>{z) + : ' 



7T 2 T 2 



f(z) h(z) 



where the primes denote derivatives with respect to z. Equation (3.9) has three (real- 
valued) regular singular points at z = — 1, l,oo and one irregular singularity at z = 
0, thus we must proceed with some care. To numerically compute the quasi-normal 
frequencies, we employ a method first developed in [31]. That is, we pull out from (p 
its infalling behavior near the horizon as well as its leading normalizable behavior near 
the boundary and define a new function <p according to 

(p(z; u>, k) = f(zy tu > /(47rT) z 5/2 e-* ka2T/z v(z; u, k). (3.10) 

The equation for <p, obtained from (3.9), does not have any singular points in the 
interval z G [1, 0]. Hence, we can approximate <p in that interval by a series expansion, 
with M terms in it, around a point inside the interval, say zq = 1/2 (so that the radius 
of convergence of the series covers the two endpoints of the interval): 

M / l\ m 

<f(z; u,k) = y~] a m {u, k) I z - ^ j . (3.11) 

m=0 ^ ' 

The equation for <p can then be cast in the following form 

M 

^2 Amn{u, k)a n (u, k) = 0, (3.12) 

m=0 

which is a matrix equation describing M + 1 linear equations for M + 1 unknowns 
a n (uj,k). For any fixed k, the quasi-normal frequencies, denoted by u n , are then ob- 
tained by solving the following polynomial equation in u: 

det A mn (u, k) = 0. (3.13) 

As we increase the size of the matrix A defined above, the poles are found to move 
around in the complex frequency plane slightly. However, depending on the desired 
precision, they converge for a large but finite value of M. For concreteness, we ran all 
our computations with M = 300, for which we found a confidence level up to fourteen 
significant digits for the imaginary part of the first mode. 

Now, since the boundary field theory under consideration is non-commutative, there 
is another scale in the problem besides the temperature, namely, the length scale \/9 
(or rather A 1//4 \/^ for our strong coupling setting). As such, the poles in the retarded 
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Green's function of the operator 0(u,k), or equivalently, the quasi-normal frequen- 
cies of the dual bulk fluctuation, depend on k, T and l/y/9 and can be conveniently 
parametrized as 

ou n (k,T,0)=Tf n ^,TV9). (3.14) 

The index n > which is sometime referred to as overtone number denotes the number 
of poles. For instance, n = is assigned to the pole with the smallest magnitude. Aside 
from special regions in the parameter space ofk/T and Ty/9, where one can determine 
the functions f n semi-analytically, f n s should, in general, be determined numerically. 
As we alluded to earlier, one expects the deviations from the commutative results to 
be more pronounced in the region of parameter space with T 3> O^ 1 / 2 and k 3> T" 1 
which also implies k 3> O^ 1 / 2 . As we demonstrate below, this expectation is confirmed 
in our numerical computations. 

The plots (a), (b) and (c) in Figure 1 show, respectively, the location of the first 
few quasi-normal frequencies for k = nT/lO, k = ttT and k = lOnT for a fixed value 
of naT = 1. The red dots show the location of the quasi-normal frequencies of ip 
on the non-extremal background (2.1) while the black dots represent the quasi-normal 
frequencies of (p on the AdS 5 Schwarzschild background. Note that even in the presence 
of a non-zero non-commutativity parameter, the poles are all located in the lower half 
of the complex frequency plane. Indeed, one can easily show that regardless of the 
range of parameters no poles will move to the upper half plane, indicating the stability 
of the system. See the appendix for a proof. 

Focusing on the dominant quasi-normal frequency (the one with the largest imagi- 
nary part), Figure 2 shows the imaginary part of u/Tasa function k/T, with T\/9 held 
fixed, for that mode for both non-commutative and commutative cases. We have also 
plotted the behavior of its imaginary part versus temperature for k = 5/100, k = 1 and 
k — 10 (in units of or 1 ) where we have fixed a = A 1 / 4 ^ = 1. Those results are shown 
in Figure 3. In each plot, the red curve shows the non- commutative result while the 
black curve shows the commutative result (9 = 0) for the same value of momentum k. 
Fitting those curves with a simple power law Imwo = —a T 7 in the regime of T\/6 ^> 1, 
we obtain the results reported in the Table I. 

From our numerical results above one can then conclude that in the regime of pa- 
rameters with both T\/6 ^> 1 and k\/9 ^> 1, the decay rate of a mode with momentum 
k into the non-commutative bath is parametrically faster compared to the decay of the 
same mode into an ordinary commutative bath. In the regime of high temperatures (as 
compared to 9~ 1//2 ), the non-local effects in non-commutative thermal baths are due 
to the large transverse size of the typical modes in the bath. These modes provide for 
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Figure 1: The red and black dots show, respectively, the location in the complex frequency plane 
of the first few quasi-normal frequencies of a minimally-coupled massless scalar on the non-extremal 
background (2.1) and the AdSs Schwarzschild background for (a) k = irT/10, (b) k — irT and (c) 
k = lOirT, where we have kept fixed iraT — 1. The plots have been generated with M — 300. 




5 10 15 

k/nT 



Figure 2: The behavior of the (negative of the) imaginary part of ui/T as a function of k/T for the 
dominant quasi-normal frequency in the commutative case (8 = 0) shown by the back curve, and in 
the non-commutative cases with iraT =1/5 and iraT = 3/5 which are shown by the red and blue 
curves, respectively. 



a rapid dissipation of the high energy modes injected into the non- commutative bath. 
In the case of a black hole or deSitter space, the event horizon appears to an observer 
who is hovering outside the horizon increasingly hotter as the observer gets closer to 
the horizon. Our analysis does show that at increasingly higher temperatures the rate 
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(c) 

Figure 3: The behavior of the (negative of the) imaginary part of the dominant quasi-normal 
frequency as a function of temperature for k = 5/100 (a), k = 1 (b) and k — 10 (c) in units of 1/V0. 
The red curves show the behavior for the non- commutative boundary theory while the black curves 
show the behavior of the same mode in the commutative theory (0 = 0). 



Commutative Non-commutative 
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a 


7 


a 


7 


0.05 


5.941 


1.318 


8.629 


1.000 


1.00 


13.398 


1.910 


8.623 


1.000 


10.00 


103.64 


1.986 


8.623 


1.000 



Table 1: Fits of the form Imw = — aT 7 for the curves shown in Figure 3 in the regime of Ty/0 > 1. 

of dissipation increases parametrically much faster in the non-commutative case. This 
appears consistent with the fast scrambling property of horizons. 

3.2 Dilaton and Axion Fluctuations 

Similar to our analysis in the previous section, we consider in this section the fluctua- 
tions of the dilaton and the axion in the background (2.1) and analyze the approach to 
thermalization of the boundary non-commutaive bath when it gets perturbed by their 
dual operators. 

To start, consider the bosonic part of the action for type IIB supergravity. In the 
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Einstein frame, ds 2 E 



\fge $//2 <is 2 tr , we can write it as 



S 




(3.15) 



where 



^3 — -^3 — CH 3 , 

F$ — F 5 — C2 A H[ 



3, 



F 3 = dC 2 , 
F 5 = dC 4 . 



(3.16) 



The type IIB equations of motion can be deduced from (3.15), supplemented by the 
self-duality condition 



In particular, the background (2.1) is a solution to the equations of motion obtained 
from the above action. 

We now proceed by considering the fluctuations of the dilaton and axion around the 
background (2.1). As we alluded to earlier, whereas there are no dual gauge invariant 
local operators in position space, it is possible to write down the linearized coupling of 
a supergravity field to the noncommutative Yang-Mills modes in momentum space as 
shown in (3.6), with O(k) being a gauge invariant operator of the type given in (3.5) 
(see, for example, [33, 34, 35]). For the dilaton and the axion, we can think of these 
operators as the non-commutative gauge invariant generalizations of Tr F 2 and Tr FF, 
respectively. 

The equations of motion that result from the above action for the dilaton $ and 
the axion C are [32] 



*F 5 = F 5 . 



(3.17) 



V 2 $ = e 2 *d u Cd»C - 




12 



12 



1 fii^p 1 1 



VV%C) 




6 



(3.18) 



where it is understood that the indices //, u, . . . run in ten dimensions. 

Now, considering small perturbations as defined by $ = $ + 4> an d C = C + x, 
where <£> and C = are the background values given in (2.1), we get the following 



differential equations in momentum space (after linearizing): 



= f{z)z 2 4>"{z) + [zf'(z) - 3f(z)]z<ft(z) + 
lQn 6 T 6 9 2 



f(z)z 2 

7T 2 T 2 



k 2 



z 6 f(z) 



[l + 2h(z)) + 



h{z) 



f{z)h'(z) Ah'(z) {2n 6 T 6 6 2 % 2 T 



2f(z)h(z) h(z) 



+ 





7T 2 T 2 


z A 


2 


z ). 







f( z y f( z ) h( z ) 

h"{z) 
f 2h{z) 



(3.19) 



and 

= f(z)z 2 X "(z) + 



znz) + Z -^±-3f(z) 



zx'(z) + 



7T 2 T 2 



UJ 



k 2 



f(z) Hz) 



X(z), 
(3.20) 



where, as before, we defined z = Uh/u and k = k^. 

Close to the horizon, the solutions to both equations behave as 



0,X ~ f(z) 



(3.21) 



where the + and — signs represent outgoing and ingoing waves, respectively. Thus, 
we discard the solutions with the + sign. Near the boundary, and for non-vanishing 
momentum, the solutions take the form 



0, X ~ A(u>, k) z * e k(>7TT/z + B(uj, k) z V fc ^ T/z , 



(3.22) 



where A = 13/2 for the dilaton and A = 1/2 for the axion. Notice that their asymptotic 
behavior is similar to that of a minimally-coupled scalar field as given in (3.8), although 
with a different characteristic exponent for z. The first term above blows up near the 
boundary so we take it to be proportional to the source for the boundary theory dual 
operator. Equations (3.19) and (3.20) have also three real- valued regular singular points 
at z — — 1, 1, oo and one irregular singularity at z = 6 . The quasi- normal frequencies 
of the dilaton and axion fluctuations in the background (2.1) can be computed using the 
same method we employed in the previous section for the case of a minimally-coupled 
massless scalar field. 

In Figure 4 we plotted the behavior of the imaginary part of the dominant quasi- 
normal frequency as a function of the temperature for k = 5/100, k = 1 and k — 10 
(in units of gT 1 ) where we fixed a = A 1//4 v^ = 1. In each plot, the blue curve shows 
the result for the dilaton fluctuation while the purple one corresponds to the axion 
fluctuation. For comparison purposes, we also plotted the commutative result (6 = 0) 
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4 6 8 10 12 

T 



(c) 

Figure 4: The behavior of the (negative of the) imaginary part of the dominant quasi-normal 
frequency as a function of temperature for k = 5/100 (a), k = 1 (b) and k — 10 (c) in units of 1/yd. 
The blue curves show the behavior for the dilaton fluctuation whereas the purple ones represent the 
axion fluctuation. For completeness we are also plotting in black the behavior of the same mode for 
the commutative case. Each point on the plot has been generated with M — 300. 



Dilaton Axion 


k 


a 


7 


a 


7 


0.05 


19.188 


1.1094 


6.9056 


1.3758 


1.00 


20.691 


1.7903 


17.047 


1.8442 


10.00 


114.31 


1.9552 


109.14 


1.9703 



Table 2: Fits of the form Imw = —aT 1 for the blue and purple curves shown in Figure 4 in the 
regime of T > 9~ 1/2 . 

in black which is same result as for the minimally-coupled massless scalar field in the 
commutative case 7 . In this case, the results for the best fits are given in the Table II. 

These results support the main conclusion found in the previous section: in the 
regimes at which non-local effects become important, i.e. T ^> O^ 1 ^ 2 and k 3> O^ 1 ^ 2 , the 
decay rate of a mode with momentum k into the non-commutative bath is significantly 

6 For k — the point z = is regular in both cases. We do not consider this case here. 
7 This is because when the background is AdS (Schwarzschid black hole), the equations of motion 
for the dilaton, the axion and a minimally-coupled massless scalar are all the same. 
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faster in comparison to the decay of the same mode into an ordinary commutative bath. 

It is interesting to note that in the regime where k 3> #~ 1//2 the behavior of the 
leading dominant pole (as a function of temperature) is almost the same for the dilaton, 
the axion, and the minimally-coupled massless scalar field, which leads to a "universal" 
thermalization time in the UV. This behavior could be attributed to the presence of the 
open Wilson lines in the definition of the gauge invariant operators in non-commutative 
gauge theories which, as a result, dominate the UV behavior of the two (and higher) 
point functions of the gauge invariant operators [27]. 

4. Discussion 

In this paper we showed that strongly coupled non-commutative gauge theories display 
at high temperature, a fast dissipation rate not seen in local field theory. At the core 
of this result is the UV-IR connection embodied in the large transverse size of high 
momenta mode. In order to assess how crucial the strong coupling ingredient is, it is 
useful to compare the behavior at high temperature in the weakly coupled regime. An 
example suggestive of this weakly coupled regime can be found in [36]. Although the 
model studied there is not a gauge theory, the leading perturbative contributions are 
similar to the ones found in gauge theories and therefore can be considered exemplary 
of the behavior at weak coupling. In that work, the authors computed the decay 
rate of a disturbance of the heat bath and found that the planar contributions to the 
thermalization rate are always larger than or equal to the contributions from the non- 
planar diagrams 8 . In Figure (5) we show an example of such behavior for k = 10 in 
l/y/6 units and small frequency co. It is clear from the figure that the enhancement of 
the thermalization rates found here in strongly coupled non-commutative baths is not 
present in the corresponding weakly coupled theories. 

It seems, then, that strong coupling is necessary to have faster thermalization 
together with the non-locality inherent of non-commutative theories. We conclude 
that this provides an important ingredient for the study of fast scramblers. Indeed, 
there have been studies of the fast scrambling conjecture in the framework of Matrix 
theory [3, 4, 5], which is known to be related to non-commutative theories [39, 40]. As 
a final remark, it would be interesting to investigate the decay rates in other non-local 
theories such as the theory dual to the near horizon limit of a stack of NS5-branes 
[28], the so-called little string theory, and compare the results with the findings in this 
paper. 

8 Planar diagrams are the same as in ordinary field theory. Non-planar diagrams have an extra 
factor that depends on the Moyal tensor e~i p ^ "' qu . For a review of non- commutative quantum field 
theories, see [37, 38]. 
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Figure 5: Comparison of the thermalization rates for the commutative (black) and noncommutative 
(red) cases at weak coupling. 
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Appendix 

In this appendix, we show that the poles in the retarded Green's function of the bound- 
ary theory operator which is dual to a minimally-coupled massless scalar field in the 
bulk, are all in the lower half of the complex frequency plane. This then implies that 
the non-commutative field theory under consideration is stable against the small per- 
turbations caused by turning on the aforementioned operator. Our analysis here follows 
the argument given in [18] where it was shown that the quasi-normal frequencies of a 
minimally-coupled massless scalar field in the Schwarzschild AdS background are all 
located in the lower half of the complex frequency plane. 

The equation of motion for a minimally-coupled massless scalar in the non-extremal 
background (2.1) is given in (3.4). Defining the tortoise coordinate w* by 



du 




1 




the equation (3.4) can be put in the form of a Schrodinger equation 







(5.3) 
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In the above expression, it is understood that u is a function of u*. Note that the 
potential blows up in the asymptotic boundary (u* — > 0) and vanishes exponentially at 
the horizon where u* — > — oo. 

Suppose now that ^ is a quasi-normal mode with the associated quasi-normal 
frequencies ui n . Since, by definition, the quasi-normal mode \l/ is infalling near the 
horizon, we isolate its near horizon behavior and define 



~ e 



(5.4) 



Note that ip(u*) vanishes in the asymptotic u* — > region. Substituting (5.4) into 
(5.2), and changing back to the w-coordinate, one easily obtains 



d u [u 2 f(u)d u ip(u)] - 2iuj n d u ip{u) +U(u)i/j(u) = 0, 



(5.5) 



with 

Now, multiplying (5.5) by ijj and integrating the result, we obtain 



(5.6) 



du 



u 2 f(u)\f(u)\ 2 + 2iu n $(u)ij'(u) + U(u)\^(u)\' 



0, 



(5.7) 



where ewe have also performed an integration by parts. Subtracting (5.7) from its 
complex conjugate yields 



i; 



du 



u n ip(u)^'(u) + u n il)(u)^'(u) = 0, 



(5.8) 



which, after an integration by parts, results in 



(Imw n ) / tp(u)ip'(u)du = --u r , 



[u 



(5.9) 



Note that, from the above equation, for Reu n ^ 0, one obtains lmu n ^ 0. Substituting 
(5.9) into (5.7), one obtains 



i: 



du 



u 2 f(u)W(u)\ 2 + U(u)\ij(u)f 



Imco> n 



0J r . 



(u = 1)| 2 = 0. 



(5.10) 



Since the potential U(u), given in (5.6), is positive definite for all (real) values of k and 
a in the range u G [1, oo), one deduces from (5.10) that Imw n < 0. 
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